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ABSTRACT 



The properties of a massive star prior to its final explosion are imprinted in the 
circumstellar medium (CSM) created by its wind and termination shock. We perform 
a detailed, comprehensive calculation of the time- variable and angle-dependent trans- 
mission spectra of an average-luminosity Gamma-Ray Burst (GRB) which explodes in 
the CSM structure produced by the collapse of a 20 Mq, rapidly rotating, Z = 0.001 
progenitor star. We study both the case in which metals are initially in the gaseous 
phase, as well as the situation in which they are heavily depleted into dust. We find 
that high-velocity lines from low-ionization states of silicon, carbon, and iron are ini- 
tially present in the spectrum only if the metals are heavily depleted into dust prior 
to the GRB explosion. However, such lines disappear on timescales of a fraction of a 
second for a burst observed on-axis, and of a few seconds for a burst seen at high- 
latitude, making their observation virtually impossible. Rest-frame lines produced in 
the termination shock are instead clearly visible in all conditions. We conclude that 
time-resolved, early-time spectroscopy is not a promising way in which the properties 
of the GRB progenitor wind can be routinely studied. Previous detections of high ve- 
locity features in GRB UV spectra must have been due either due to a superposition 
of a physically unrelated absorber or to a progenitor star with very unusual properties. 



1 INTRODUCTION 



For a few hours after their onset, the afterglows of Gamma- 
Ray Bursts (GRBs) are the brightest sources in the far Uni- 
verse. Their high luminosity, together with their powerlaw, 
featureless spectrum, make them ideal sources to probe their 
surrounding environment through the absorption lines im- 
printed in their spectra. 

Theoretical studies have suggested that long Gamma 
Ray Bursts (GRBs) are produced by the collapse of rapidly- 
rotating, chemically homogeneous, massive stars (e.g. Mac- 
fadyen & Woosley 1999). The association between GRBs 
and massive stars, which has been observationally supported 
(Stanek et al. 2003 Hjorth et. al||2003 ), makes absorption 



et al. (2008); see also Whalen et al. 2008 for an extended dis- 
cussion on absorbers in GRBs), yielding constraints on the 
nature of the absorbing medium. An especially well studied 
burst was GRB 021004, whose high-resolution spectroscopy 
revealed a complex velocity structure of the absorbing ma- 
terial, with velocities up to > 3000 km/s. The interpreta- 
tion of these features has been controversial. Initial studies 
claimed that the high velocity lines were a direct proof of the 



association of GRBs with WR stars. |Starling et al. (2005) 



argued that the lines must be produced in a fossil stellar 
wind with hydrogen enrichment from a companion. |Mirabal| 
et al. ( 2003 ) and Schaefer et al. ( 2003 ) interpreted the lines 



studies potentially useful as a new way to probe star-forming 
regions at intermediate and high redshifts, and/or the last 
hundreds of years of the progenitor evolution. Absorption 
lines imprinted by the material ejected by the star prior to 
its explosion allow one to probe the velocity structure and 
metal content of the eject a. High-resolution spectroscopy of 
GRB afterglows has been performed in a number of cases 



as the result of shells of material which are present around 
the progenitor. 

More recent studies have however cast doubts on the 

([2006) performed a de- 



initial interpretation. Lazzati et al 



(e.g. Moller et al. (2002); Matheson et al. (2003); Mirabal 


et al. (2003); Schaefer et al. (2003); Starling et al. 
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2005 


); 


Fiore et al. (2005); 


D'Elia et al. 
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2009 


); Chen et al. 
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2008 


); 



Prochanskaet al. (2008a), (2008b); Fox et al. (2008); Thoene 



tailed time dependent analysis, taking into account the burst 
flash ionization. Even though they still considered the wind 
of the WR progenitor star as the best absorber candidate, 
they pointed out that such an interpretation would require 
a termination shock at a distance of at least 100 pc. Such a 
large radius of the termination shock is somewhat at odds 
with current wind models and could be accounted for only 
if the progenitor were an extremely massive star evolving in 
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a fairly low-density environment. Chen et al. ( 2007[ ) stud- 
ied a sample of 5 GRBs with high velocity features (among 
which GRB 021004). Using fine-structure transitions and 
the presence of low-ionization species, they argued that the 
location of the absorber is very likely at a large distance 
from the burster (> 1 kpc), favoring an intervening halo, 
physically unrelated to the GRB, as the location of the ab- 
sorbing material. Similarly, Vreeswijk et al. (2007) modeled 
fine structure lines present in the spectrum of GRB 060418 
with a UV pumping model and concluded that the absorbers 
are at a distance of 1.7 ± 2 kpc. These results are consis- 
tent with a GRB completely ionizing the local GSM, and 
with high- velocity absorbers elsewhere in the host galaxy or 
intervening intergalactic space (see also Prochanska et al. 
(2008a,b) and D'Elia et al. (2009a) for similar findings in 
the case of other GRBs they studied). 

Most of these studies do not take into account the pos- 
sible shielding effects of an unusually dense GSM, as well as 
any effects due to dust. Additionally, when photoionization 
calculations are performed, most assume a steady state so- 
lution, whereas the radiation fields of GRBs are highly tem- 
porally variable, with initial burst timescales of a few to tens 
of seconds, and the afterglow varying over minutes to hours. 
Disentangling the physical origin of the spectral absorption 
features seen in GRB afterglows is clearly of great impor- 
tance for a full understanding of the GRB phenomenon. 

The aim of this paper is a detailed, systematic study 
of the time-dependent transmittance spectrum of a GRB 
exploding in a dusty circumburst medium (GSM) shaped 
by its progenitor star. Recent studies have suggested that 
long GRBs are produced by the collapse of rapidly-rotating, 



chemically homogeneous, massive stars (Yoon & Langer 
2005) IWoosley Heger|[2006l |Yoon et al.||2006| |CantieIlQ 



et al. 120 07*). This kind of evolution is believed to take place 



at low metallicity, where stellar winds are weak and angular 
momentum loss by winds is marginal. The wind of numerous 
types of GRB and SN progenitors has been simulated by |van| 
Marie et al. (2008), including one model which recreates an 



anisotropic wind close to the star. This theoretical scenario 
has found support in recent observations. Gamp ana et al.| 
( |2008 ) performed a detailed study of the X-ray spectrum of 
GRB 060218, discovering a larger than normal N/0 ratio 
in the surrounding of the burst progenitor. They concluded 
that only a progenitor star characterized by a fast stellar 
rotation and sub-solar initial metallicity could produce such 
a metal enrichment. 

In this paper we perform a detailed and comprehensive 
study of the effects of the GRB and afterglow radiation on 
the GSM structure produced by the stellar wind of a rapidly- 



rotating, low- metallicity, massive progenitor star (van Marie 
et al.|[2008 ). We aim at predicting the time- variable, angle- 



of-sight dependent transmittance spectra of a GRB after- 
glow which resulted from the collapse of such a star. We 
study two cases, that of a GSM whose metals are initially in 
the gaseous phase (i.e. no dust), and that of a GSM whose 
metals are heavily depleted into dust prior to the GRB ex- 
plosion. Even though the formation and survival of dust in 
the wind of Wolf-Rayet stars may seem surprising, dust is 



observed in WR winds (Allen et al. 1972 Williams et al 



1990). The gas temperature and the exposure to the star's 



the gas phase, the grains would quickly sublimate. One key 
ingredient that allows the dust to survive is its capability 
of radiatively cooling to temperatures lower than the gas in 
which it is embedded (e.g. (Lazzati 2008)). Even so, there 
are still many unknown riddles in the way in which dust par- 
ticles are born and survive in WR winds, but their presence 
has been clearly detected (e.g. Growther 2003). 

Although it is presently unknown whether the partic- 
ular stellar evolution model that we consider here is truly 
representative of critically rotating stars in general, we be- 
lieve that the ability to compare theoretical absorption fea- 
tures with observational data can be potentially useful to 
the study of both GRBs and massive star evolution. 



2 THE CSM OF THE PROGENITOR STAR 
AND SPECTRA CALCULATIONS 

2.1 The CSM of the Progenitor Star 

Van Marie et al. (2008) computed the evolution of the cir- 
cumstellar density profile around a rapidly rotating, 20Mq, 
low metallicity star (Z = 0.001) star, which is considered 
representative of a typical GRB progenitor. Due to the 
anisotropy associated with the rotation of the star, the cal- 
culation yielded density profiles which are dependent on the 
angle between the rotation axis (presumably coinciding 
with the GRB jet axis), and the line of sight. 

We choose 3 different radial profiles to show the angu- 
lar and radial dependence in detail, at = 0, 25, and 80 
degrees. The calculations of [van Marie et al.| (|2008) provide 
the density, the radial and tangential velocity Vr and ve , the 
ionization state of hydrogen, and the internal energy as a 
function of radius r and angle 0. The internal energy is used 
to compute the initial temperature as a function of radius by 
means of the ideal gas equation of state with the polytropic 
index 7 = 5/3. The initial ionization state of elements heav- 
ier than hydrogen is computed using the photoionization 
code of |Perna &: La zzati' ('2002), assuming photoionization 
equilibrium under the blackbody radiation field from a star 
with a characteristic temperature of 1.4 x 10^ K and a lu- 
minosity of 4 X IO^Lq, which correspond to the properties 
of the star at the end of its life as given by |van Marie et al.] 
( [2008 p . 

The density profile is assumed constant in time through- 
out the calculation, which is reasonable since the hydro- 
dynamical timescales are much longer than the ionization 
timescales. For the same reason, the radial velocities of the 
gas are assumed constant in time as well. 

For the metal content of the GSM, we use the metallic- 
ity history of the stellar photosphere as calculated by | Yoon] 
et al. (20061). This, combined with the radial velocity infor- 
mation, is used to compute a radially-dependent metallicity 



UV radiation make the wind a very harsh environment for 
dust particles. If the dust were to share the temperature of 



profile. The outer shell, with near-zero radial velocity, has 
the initial metallicity (75% H 25% He with trace elements), 
while the inner region, dominated by the stellar wind, is 
very helium enriched. The helium enrichment is due to the 
rapid rotation of the star; it causes the star to evolve homo- 
geneously because of mixing, so that at the end of the star's 
life nearly all of the hydrogen has burned into helium. 

For the species that are not included in the calculations 



of Yoon et al. (2006) (Sulfur, Argon, Galcium, and Nickel), 
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Table 1. Abundances used for the outer shell and the inner re- 
gion. These are given on a base 10 logarithmic scale normalized 
so that the abundance of hydrogen is 12. 



Element Outer Shell Abundances Inner Region Abundances 



He 10.90 13.84 

C 7.28 8.64 

N 6.68 10.33 

O 7.61 8.62 

Ne 6.75 8.92 

Mg 6.17 8.47 

Si 6.20 8.52 

S 5.9 5.9 

Ar 5.6 5.6 

Ca 5.0 5.0 

Fe 6.19 8.51 

Ni 5.0 5.0 




Figure 1. Temperature, Radial Velocity, and Density Profiles. 
The solid, dotted, and dashed lines represent the ^ = 0, 25, and 
80° lines of sight, respectively. 



we assumed a metallicity of Z — 0.001 with the relative 
weights of the species as in the solar case ( [Anders Grevesse 
|1989 ), for both regions of the CSM. The abundance values 
we used for all 12 species are given in Table[l] The HI column 
density is about 1.5 x 10^^ cm~^. 

The velocity, density and temperature profiles produced 
by the simulations above are displayed in Fig[l]for the three 
lines of sight under consideration. Full 2D plots of the veloc- 
ity, density, and temperature are shown for completeness in 
Fig [2] [3] and|4] These figures display the physical conditions 
of the CSM just as the star is about to explode. 

An component which is not included in these models is 
dust. As a matter of fact, while it is known that dust is pro- 
duced by stars, there is very little information on the precise 
radial distribution, grain-size distribution and composition 
of the dust grains in the close CSM of the star. Therefore, in 
our study we consider two limiting situations, one in which 
there is no dust (all the metals are in the gaseous phase 
to start with), and one in which the metals are heavily de- 
pleted into dust at all radii. In this case, we take a dust 
population similar to that of the Milky Way, consisting of 
graphite and silicate grains with a power-law distribution 
in size with slope -3.5, and a range of sizes between .005 
jim. and .25 jim.. Since the metallicity of the CSM is radi- 
ally dependent, we determine at each radius the dust-to-gas 
ratio by assuming that 90% of the least abundant element 
among the ones making up the grains (C, Mg, O, Si, Fe) 
is depleted into dust. We point out that the precise details 
of the initial dust distribution are not very important for 
our results, since the dust grains close to the GRB loca- 
tion are quickly destroyed by the intense radiation from the 
burst. Similarly, any uncertainty in the destruction rates by 
a factor of a few would not affect sensibly our results, since 
recycling is extremely fast with respect to photoionization. 
However, there is a fundamental difference between an initial 
condition with no dust and an initial condition with dust. 
In the former, low-ionization species such as CII, Fell, Sill 
are photoionized by the strong radiation field of the pro- 
genitor star in its close CSM, where the high- velocity, free- 
streaming wind resides. Therefore, high-velocity lines from 
these elements are never visible. In the latter case, these el- 
ements are initially shielded within the dust grains. Then, 
when the star explodes and the prompt GRB and afterglow 




Figure 2. Radial Velocity Map: The radial velocity map for the 
polar region between angles ^ = and 6 = rad. The color 
scale is in units of log [km/s]. The 0-10 pc region has velocities 
ranging between 2730 and 5332 km/s, with the highest speeds 
found at smaller polar angles. The over-plotted white dash lines 
indicate the lines of sight displayed in Fig^ Note that one line of 
sight lies along the pole (y axis) and that positive velocity here 
indicates outfiow. 



radiation impinge on the grains, those depleted elements are 
recycled into the medium, and they are visible for a period 
of time until they get photoionized. Thanks to the use of the 
code by Perna & Lazzati (2002), we are able to model self- 
consistently the destruction of the dust grains, their redistri- 
bution into the CSM, and their subsequent photoionization. 
When there is dust, high- velocity lines from low-ionization 
species become visibl^ 

We thus have all the initial conditions of the medium at 



^ We need to point out that the results we display in the following 
would change qualitatively if the dust composition were radically 
different from the MW distribution assumed here. For example, if 
the grains, instead of being composed of Fe, Mg, Si, O and C (as 
in the MW dust distribution), were made up of other elements, 
then lines from those elements would be seen coming from the 
wind. 
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1 



The simulation of time-dependent photoionization and 
radiation transfer was performed using the cod e by [Perna 



l-?s 



1-21 



Figure 3. Density Map: Similar to Fig. |2] but for the density of 
the medium. The color scale is in units of log [g/cm*^]. 




Figure 4. Temperature Map: Same as in Figs. |2] and [s] but for 

the temperature of the CSM prior to the star explosion. The color 
scale is in units of log [K] . 



the time at which the star explodes and produces a GRB: 
temperature, density, radial velocity, metallicity and dust 
as a function of radius and polar angle. These conditions 
provide the starting point for the computation of the time- 
dependent transmittance spectra expected during the GRB 
afterglow phase. 



2.2 Numerical Setup and Radiative Transfer 

We divided each line of sight (the three chosen values of 
0) into 1909 radial bins, matching the resolution of the data 
generated by [van Marie et ah ([2008). In the time coordinate, 
we chose logarithmically spaced time steps starting at 10~^ 
seconds and ending at a time of 3 days. It's worth noting 
that time t = is defined as when the first photons reach 
the observer. Thus, observational effects like superluminal 
velocities of ionization fronts are apparent, not actual. Ion- 
ization profiles are, instead, in the rest frame of the burst. 
Complete ionization profiles (necessary for generating the 
spectra) were therefore extracted as a function of time. 



& Lazzati' ("2002) (see also Perna et al. (2003); Lazzati & 



Perna (2003)). This code keeps track of the 13 most abun- 
dant elements (H, He, C, N, O, Ne, Mg, Si, S, Ar, Ca, 
Fe, Ni) and computes their ionization state as a function 
of time. The continuum opacity is also part of the compu- 
tation, while the opacity due to resonant lines is calculated 



separately (see ^2.4). The radiative transfer is performed in 
the 1 eV — lOkeV spectral range. 

We note that, once the ions from the sublimated grains 
return to the gaseous phase, they are given the temperature 
of the gas at that location. This is likely to be inaccurate, 
but unfortunately, at this stage it is not possible to imple- 
ment a more accurate prescription in our code for both a 
technical and a theoretical limitation. Technically, once the 
ions are ejected from the dust, they are not distinguishable 
(in the code) from the gas ions. The code does not allow for 
two components of the same ion at different temperatures. 
Furthermore, even if the code were to include this possibility, 
there is not a single dust temperature. The temperature of 
the sublimating dust grains depends on the grain size and on 
the distance to the burst. Therefore one would need to con- 
sider a continuum of temperatures and the problem would 
become numerically untreatable. Finally, there is a theoreti- 
cal uncertainty in what is the distribution of the ejected ion 
energies. One possibility is that the ions are ejected with a 
thermal energy distribution corresponding to the tempera- 
ture of the dust. The issue with this procedure is that the 
ejected particles have to overcome a potential barrier (the 
binding energy to the grain) to leave the dust grain, and 
there is no clear recipe on how to compute the ejected par- 
ticle energy distribution. For big grains, it is likely to be an 
exponential tail (and therefore non-thermal), but for small 
grains quantum effects are important. We consider a de- 
tailed implementation of all these effects, albeit interesting, 
beyond the scope of this paper. 



2.3 The Light Curve of the GRB and its 
Afterglow 

The light curve of the GRB and its afterglow was modeled 
with the sum of two separate components: a prompt emis- 
sion component and an afterglow component. The models 
we adopted are described in this section. For each, we let 
our light curve take the form: 



(t, V, 0) 



AT{t)F.{u) (^^y 



(1) 



where A is an appropriate normalization constant, T(t) con- 
tains the time dependence of the radiation, Fu{iy) is the flux 
at a given frequency and represents the half opening 
angle of the jet core within the context of a structured jet 
dRossi et al.|2002| [Zhang Meszarosf2002[ ). This is observed 
to vary from less than 1° to as much as 15° (Panaitescu 2007; 
Liang et al. 2008) and is taken here to be = 2.5° . Note 
that the angular dependence of the luminosity does play a 
significant role in the ionization timescale of the various el- 
ements at different polar angles. 

Both the spectral function Fjj(y) and the time function 
T{t) are normalized so that the integrals over frequency and 
time, respectively, are equal to 1, making the normalization 
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of the luminosity entirely contained within the constant A. 
Additionally, we note that 



(2) 



where Lis© is the effective isotropic luminosity for the line 
of sight we are interested in, an important input onto the 
radiative transfer code. 



2.3.1 The Prompt Emission 

For the GRB prompt emission, we choose a simple Gaussian 
function with a FWHM of 10 seconds and a peak emission 
at 10 seconds ( Kouveliotou et aL]|1993 ): 



with P = 0.975 and Aagi a fitting constant, picked so that 
at peak emission Tag ^^peak) ~ ^' "^^^ powerlaw break 
is at tb = 100 s. We analytically find the peak time to be 
^peak ~ 132.677 s, and our fitting constant was determined 
to be AAGt = 1.72796. 

To normalize the frequency dependent component we 
use the observation that R — 15 for a z = 1 GRB. (van 



(3) 



Tp{t) = AptG[t, = 10 s, cr = 10/(2V21n2) s] . 

The value of the constant Apt — 1.00935 s~^ is determined 
so that the Gaussian is normalized to 1 in the interval < 
t < oo. 

We use the broken power law form given by [Band et] 
al. (1993) for the frequency dependent component of the 



[Paradijs et al.|2000| ). To fit this, we used the R-Band Bessel 
filter zero point: = at a flux density of 3080 Jy at 
= 4.69 X 10^^ Hz. This gives an = 15 flux den- 
sity of .00308 Jy at the same frequency. Using a luminosity 
distance dL{z = 1) = 6626.5 Mpc (for a cosmology with 
Qa 0.73, Qm 0.27, and Hq 71 km/s/Mpc), we find 
a specific luminosity ^z^(^peak' ^r) ~ ^•^'^ ^ lO^^erg/s/Hz. 
Since TAG(tpgak) = 1' ^^en i^^(tpeak' ^r) = ^^(^r) = 
AAGfT^^ = 1.62 X 10^^ erg/ s/ Hz with the typical value of 
a = 0.65 (corresponding to an index p = 2 of the electron 
distribution). This readily yields AAGf = 5.57 x 10^^ Thus, 
the spectral component is 



prompt emission: 



FAGiiy) = 5.57 X 10^V"°-^^ erg s" 



^Hz- 



(7) 



giving the following model for the afterglow component 

2 



(8) 



FEiE) = Apf 



{a - l3)Eo 
lOOkeV 



Q-/3 



exp (/3 - 



V lOOke V 



)' 



{a-l3)Eo ^ E 



(4) 



where E — hv \s in units of keV, [a — I3)Eq is the knee of 
the power law, taken to be 300 keV, a and /3 are the slopes 
of the power law for E < Eq and E > Eq, respectively, and 
Apf is a normalization constant. In our model, a = and 
/3 = — 2. Using these values makes the integral easy to solve 
analytically: 



poo 
J2E0 



exp(-^/300 key)dE+ 



Viooy 



with Fag{t^) and TAG{t) defined as we just discussed. This 
corresponds to a total effective isotropic energy output 
across the range 1 eV to 10 keV of about 10^^ ergs, about 
10% of the prompt emission. 



2.4 Computation of the Transmittance 

During the radiative transfer calculation, we keep track of 
the time-dependent density of each ion present in each ra- 
dial bin, hence obtaining the evolving ionization profile as 
a function of radius for each ion of interest. We track the 
following ions: HI, Hel, Hell, CII, CHI, CIV, Sill, SilV, 
Mgl, Mgll, Fell, and OH. This, combined with our radial 
velocity and temperature profile, is sufficient to generate a 
time-dependent transmittance, which we compute at infinite 
instrument resolution. 

To calculate the transmittance, we compute the Voigt 
line profile function for each individual line. 



yielding the normalization constant Apf — 2.935990 x 
10"^ keV"\ 

Thus the model for the prompt emission is 

L.(t,^,e) = £;i3oF^(^)T(t)(^)' (5) 

where E^qq{0 / Oq)'^ is the total isotropic energy output along 
the line of sight determined by 0. It is related to L^qq{0) 
through ^jgQ = ^iso(^' ^o)dt. 

2.3.2 Afterglow Model 

The time component of the afterglow emission is modeled 
as a smoothly broken power law 

Aagi 



H(a, u) 



where H{a,u) is the Voigt function: 



/o 
-c 



dy 



a"^ -\- {u — yY 



(9) 



(10) 



with a — V/A.ti/\vd and u = — uq) / Aud- Here Aud is 
the Doppler width and F is the transition rate, which is 
approximated by 7, the spontaneous decay rate jRybicki & 



Lightman|1979 ). The cross section for a given line is 



-/i; 



(11) 



TAG{t) 



(6) 



from which we can calculate the optical depth r{u). The 
transmittance is thus 



Tr{iy) = Tro{iy) * e"^. 



(12) 
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where Tro is its value before being affected by the opacity 
of the given hne. 

The transmittance was calculated at each radial bin tak- 
ing into account 21 resonant transitions from the above men- 
tioned ions; thanks to the inclusion of the Doppler effect, any 
structure present in the velocity profiles is correspondingly 
reflected in the simulated spectra. We limited ourselves to 
21 resonant transitions since this was the most computa- 
tionally expensive part of our study. The lines were selected 
based on observed GRB spectra ( |Fiore et al. 2005[) rest wave- 
length, and oscillator strength /12. The lines present in our 
calculations are listed in Table |2] The reported /12 values 
were calculated from the Einstein A coefficients given by 



Ralchenko et al. (2008). It is worth noting that we only 



treat lines that arise from resonant transitions in the ground 
state of each ion, although fine-line transitions are com- 
monly observed in GRB afterglow spectra (Vreeswick et al. 
2003, 2007; Prochanska et al. 2006; D'Eha et al. 2009a,b). At 
present time, there is no time-dependent, radiative transfer 
code that treats self-consistently transitions from both the 
ground and the excited levels, as well as the time-dependent 
evolution of the dust grains. However, we estimated the ef- 
fect of our approximation by using a time-dependent code 
(van Adelsberg et al. in prep.) which computes the evolu- 
tion of the excited states of the ions under the influence of a 
time- variable radiation field. By running this code with the 
luminosity function in §2.3 and viewing angle ^ = 0°, we 
found that, for ions such as Sill and Fell, in the innermost 
region of the wind 1 pc) the ground level depopulates 
to about 20% on a timescale of several seconds, while in 
the outermost regions of the wind 10 pc), it depopulates 
to about 30% on a timescale of several hundred seconds. 
The depopulation timescales become longer at larger view- 
ing angles. Therefore, given the timescales over which the 
transient high- velocity features live (§3), the correction to 
the strength of the ground-level resonant lines due to the 
depopulation of the ground level is not expected to be sig- 
nificant. In the case of CIV, we found the population of the 
excited levels to be negligible across the full extent of the 
wind. 



3 RESULTS 

Our calculations were performed both for a non-dusty CSM 
and for a dusty CSM. In the non-dusty case, most of the 
inner region is completely or near-completely ionized be- 
fore the burst occurs due to the ionizing photons from the 
progenitor star, vastly reducing the signatures of the high- 
velocity wind close to the star. These features, that is high- 
velocity lines from low-ionization elements, are instead visi- 
ble in the case in which metals are initially locked into dust 
grains and later recycled into the ISM by the GRB radi- 
ation. The absorption features from the outer shell, on the 
other hand, are very similar for the non-dusty and the dusty 
scenarios. Therefore, we present in the following detailed re- 
sults only for the case of a GRB exploding in the dusty CSM 
of a massive star. 

Figures|5]|6] and|7|display the results of the time depen- 
dent ionization profiles for HI, Hell, Sill, Silll, SilV, CII, 
CIV, and Fell, along lines of sights corresponding to ^ = 0°, 
25°, and 80°. To incorporate velocity information, we show 



a shaded column density map for each line of sight and a few 
key ions, with color/shade indicating the observable column 
density along the line of sight, the y-axis indicating the ve- 
locity, and the x- axis indicating time. A negative velocity 
here represents a blue-shifted velocity towards the observer. 
Note that the maximum velocity does vary with angle, as 
a result of the anisotropic nature of the CSM produced by 
a rotating star. The displayed times range logarithmically 
from about 1 x 10~^ seconds to 2.5 x 10^ seconds. 

There are some common features that are interesting to 
note. In particular, the column densities of the gaseous com- 
ponents of the low-ionization states of C, Fe and Si, which 
are locked into dust grains before the burst goes off, are neg- 
ligible at t = (explosion time), they then increase as the 
grains are destroyed by the burst radiation and recycled into 
the ISM, and eventually decrease as the ions are photoion- 
ized. This time evolution sequence is much faster at small 
viewing angles, given the stronger radiation field from the 
GRB. 

Figures [s] through 10 show several snapshots (between 



t ^ 10~^ s and t ^ 10 s) of the normalized transmission 
spectra for each of the three representative lines of sight. 
The times in these figures correspond to the vertical dotted 
lines in the aforementioned column density plots. We do not 
include information at later times because the spectra have 
stopped evolving at this point, as is apparent from Figures 
|5l[S] and|7l 

Additionally, we include cross sections across the col- 
umn density plots for the time steps shown in the sample 



spectra (Figures 11 12 and 13). These simply show, perhaps 



in a simpler manner than the column density maps, the col- 
umn densities vs. the velocity of the gas for a few selected 
ions. We indicate the most notable features in section EH] 



3.1 Manifestation of the Velocity Structure. 

To keep these results in context, it is helpful to refer to Figs. 
[l]and[2l since the velocity structure depends on 6, and also it 
is helpful to note the location of certain velocity structures 
since this has a direct impact on the ionization timescale. 
For example, along all lines of sight, the free streaming re- 
gion extends from our minimum radius of 3 x 10^^ cm, and 
experiences a slow variation in velocity from there to a dis- 
tance of about 5 X 10^^ cm. This is the innermost area of the 
CSM in our simulation. The inner subset of this region, up 
to about 1 X 10^^ cm, has the highest number density, and a 
temperature around 5 x 10^ K (Figure [ij. This is the region 
where the high-velocity features seen in the spectra are pro- 
duced. For example, for ^ = 0° (Figure [sj, this accounts for 
the 3500 km/s feature for HI, as well as the ~ 5000 km/s 
component present in several of the other ions. 

As Figure [1] indicates, there are indeed isolated struc- 
tures at specific velocities. The most dominant one is per- 
haps the least interesting, the outer shell which is nearly 
at rest {vr < 25 km/s), and is expected by the most ba- 
sic stellar wind models. This is observable as the high den- 
sity feature at the bottom of each shaded column density 
plot. Equally ubiquitous is a large bubble of shocked wind, 
which is visible in Figure [l] as the sudden drop in velocity 
at about r = 1.5 x 10^^ cm. This region occupies the ma- 
jority of our simulation space, and because of its very hot. 
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Table 2. Lines and Oscillator Strengths of the 21 Transitions considered for the calculation of our transmittance spectra. 



Ion 


Ao[A] 


/12 


Ion 


Ao[A] 


/12 


HI 


1215.67 


0.278 


CII 


903.962 


0.336 


CII 


1334.53 


0.127 


cm 


977.03 


0.759 


CIV 


1548.20 


.190 


CIV 


1550.77 


0.0952 


Mgl 


2852.50 


1.80 


Mgll 


2796.35 


0.566 


Mgll 


2803.53 


0.303 


Sill 


1190.42 


0.139 


Sill 


1193.29 


0.574 


Sill 


1260.42 


1.224 


Sill 


1304.37 


0.0928 


Sill 


1808.01 


0.00248 


Silll 


1206.50 


1.67 


SilV 


1393.76 


0.512 


Fell 


2344.21 


0.11 


Fell 


2374.46 


0.028 


Fell 


2382.77 


0.38 


Fell 


2586.65 


0.065 


Fell 


2600.17 


0.22 









^ 1 X 10^ — 1 X 10^ K temperatures, none of the ions that 
we track are present in this region. 

We present the spectra for four different ions where we 
seek to determine if the velocity structure becomes visible. 
First, we note that Hydrogen is highly ionized in the in- 
ner region and thus is a poor tracer of velocity structure. 
However, the outer shell does contribute a large hydrogen 
component and the near 1215.67 A rest frame Ly-a feature 
can thus act as a diagnostic of the outer shell. 

In the ^ = case, (Fig. |8|, we see some higher velocity 
components for ions associated with dust such as Silll, Fell, 
and CIV. These appear only at early times, when the burst 
initially destroys the dust, releasing Si and Fe. However, 
at this head-on angle, the intensity of ionizing radiation is 
high, quickly ionizing the released ions and destroying this 
signature of the high velocity wind almost immediately. 

In the — 25° line of sight, the rest frame Ly-a line 
exhibits a behavior similar to the ^ = 0° case, but on longer 
timescales. Rest frame and near rest frame features also ex- 
hibit a very similar behavior as in the ^ = 0° case. Of partic- 
ular note is Figure [To] corresponding to the — 80° viewing 
angle, in which we see finely structured Fell and Silll lines 
due to the release of gas from dust in the high velocity re- 
gion of the CSM. These survive on the order of a few seconds 
before being ionized to higher states. 



4 SUMMARY AND DISCUSSION 

In this paper we have presented a comprehensive calcula- 
tion of the spectra of a typical GRB which explodes in the 
CSM created by a massive, chemically homogeneous, low- 
metallicity rotating star. We have explored both a dusty 
CSM as well as a dust-free one. These stars are strong candi- 
dates for the typical progenitors of long GRBs; therefore, the 
models presented here, in combination with an average GRB 
luminosity and spectrum, are expected to yield the typical 
expectations for the spectra during both the prompt and 
the afterglow phase of the GRB. We emphasize, however, 
that an accurate prediction can be obtained only from a tai- 
lored time-dependent simulation for each particular burst, 
and hence our results should rather be seen as an example 
of the overall evolution, whose details will then depend on 
the spectrum and luminosity of the specific burst. 

We have performed the calculation of the transmittance 
for different lines of sights with respect to the GRB jet axis 
(assumed to coincide with the rotation axis of the star). 



thus exploring the latitudinal variation of the spectra, which 
include the ^-dependence of both the CSM profiles as well 
as of the GRB luminosity. 

The highest velocity components of the CSM, reaching 
velocities above 5000 km/s, exist predominantly at small 
viewing angles (with respect to the rotation axis of the pro- 
genitor). This is where the burst is most luminous, and thus 
ionization timescales are very fast. We do see the presence 
of several ions at these velocities, but this region becomes 
completely ionized in a fraction of a second, and an actual 
observation of the high velocity features seems hopeless, at 
least for the star considered in this study. 

Lines of sight at higher latitude are characterized 
by smaller maximum velocities, with the maximum being 
around 2800 km/s in the case of a line of sight at 80°. In our 
model, even at these latitudes, ions with lines around 1000 
km/s cannot survive the entire burst without being pho- 
toionized. From a purely theoretical prospective, the time 
variability of the high velocity lines, and in particular, the 
transition from a nearly featureless spectrum to high column 
density lines back to a featureless or weak lined system like 
the one observed in our study, is a clear signature that the 
lines are associated with the CSM, and furthermore, that 
they are due to the destruction of dust in the CSM created 
by the progenitor star. Practically, however, observational 
detection of such variability with current instruments is not 
feasible. 

Observationally, high-velocity lines have been reported 
in GRB 021004 by Fiore et al.|(|2005| after about 16 hr from 



the trigger. From the analysis presented in this paper, we 
find that these features, if associated with the free- streaming 
wind of a progenitor star like the one discussed here, could 
have survived till these late times only if the burst had a very 
low luminosity. Since this is not supported by observations 
(the luminosity of GRB 021004 was rather high), the only 
way to associate these lines to a wind is by requiring a very 
special event, in which the progenitor star would possess an 
unusually strong wind, and it were to explode in a very low 
density medium. Alternatively, the observed late-time, high- 
velocity lines could be simply due to intervening absorbers 
(Chen et al. 2007). 

In summary, during the initial moments of the GRB 
explosion, the optical absorption spectra present a very rich 
phenomenology due to the different speeds at which the wind 



components move away from the progenitor side ( van Marie 



|et al.||2008| . If we were able to observe these early phases. 
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we could derive a wealth of information on the properties of 
the progenitor stars of GRBs. More generally, this would be 
of great importance to the study of the final stages of stellar 
evolution. Observationally, however, it is extremely challeng- 
ing to catch variability on such small timescales, shorter even 
than the time it usually takes to localize the burst. There- 
fore we conclude that, if the massive stars progenitors of 
GRBs are indeed of the kind studied here (as suggested by 
theoretical investigations), then the high- velocity lines pro- 
duced in their winds are practically invisible, even for very 
low- luminosity bursts (such as the = 80° case that we 
considered). The only observational signatures of the GRB 
progenitor stars are imprinted in non-variable absorption 
lines produced in the termination shock of the wind. Previ- 
ous detections of high- velocity lines must therefore be due to 
intervening absorbers, unless the GRB progenitor star had 
extreme properties well outside of the range considered here. 
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Figure 5. Observable column density at different velocities versus time at ^ = 0°. Each panel corresponds to a different ion. They are, 
in order from top to bottom, HI, Hell, SiH, SiHI, SilV, CH, CIV, Fell. The color/shade indicates the observable column density. The 
y axis indicates at what velocity the lines will be observed. Each horizontal line element corresponds to a velocity interval of 54 km/s. 
The X axis shows time. The vertical dotted white lines indicate where we show sample spectra. 




Figure 6. Similar to figure [s] but for 6 = 25°, and with each horizontal line element corresponding to a velocity interval of 44 km/s. 
Again, the different ions are, from top to bottom, HI, Hell, Sill, Silll, SilV, CII, CIV, Fell. 




Figure 7. Similar to figurejSjbut for 6 = 80°, and with each horizontal line element corresponding to a velocity interval of 28 km/s. The 
different ions are, from top to bottom, HI, Hell, Sill, SiHI, SilV, CII, CIV, Fell. 
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Figure 8. Transmittance versus velocity (100 km/s) for Fe II 2382.04 A (left), Silll 1206.50 A (left-middle), CIV 1550.77 A (right- 
middle) and HI 1215.67 A (right) at ^ = 0° and at times t = 10"^, 0.17, 1.80, 20, and 1040 seconds (from top to bottom). Each spectrum 
is normalized and is scaled identically. 
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Figure 9. Similar to Figure [s] but for 9 = 25° 
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Figure 10. Similar to Figure |9] but for 8 = 80°. 
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Figure 11. Number Density [log cm ^] versus velocity [km/s]. Each panel represents a cross-section of Figure [s] at the times t = 0.17 s 
(solid), 1.80 s (dots), 20 s (dashes), and 21000 s (dashes-dots). The ions presented are HI (lower left). Si II (lower right), SilV (middle 
left), CII (middle right), CIV(upper left), and Fell (upper right). 
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Figure 12. Cross sections of Figure [6] of the same ions and at the same times as in Fig. |11| 
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Figure 13. Cross sections of Figure [7|of the same ions and at the same times as in Fig. 



